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Abstract 

The recent results from T2K and MINOS experiments point towards a relatively large 
value of the reactor angle #13 in the lepton sector. In this paper we show how a large 
6>i3 can arise from the charged lepton sector alone in the context of an SU(5) GUT. In 
such a scenario (tri-)bimaximal mixing in the neutrino sector is still a viable possibility. 
We also analyse the general implications of the considered scenario for the searches of 
CP violation in neutrino oscillations. 



1 Introduction 

Recently the T2K collaboration reported pQ evidence at 2.5a for a non-zero value of 
the reactor angle #13 in the neutrino mixing matrix. Quantitatively, it was found that 
0.03(0.04) < sin 2 29 13 < 0.28(0.34) at 90% C.L. for Am| 2 = ( +2.4x 10" 3 eV 2 , sin 2 26> 23 = 1 
and 5 = 0, Am| 2 and #23 being the atmospheric neutrino mass squared difference and mix- 
ing angle, and 5 being the Dirac CP violating phase in the neutrino mixing matrix (see, e.g. 
[2]). Under the same conditions the best fit value of #13 obtained from the data analysis 
is relatively large: sin 2 2^13 = 0.11 (0.14). Subsequently the MINOS collaboration also 
reported evidence for a relatively large value of #13, although with a smaller statistical 
significance [3]. A global analysis of the neutrino oscillation data, including the data from 
the T2K and MINOS experiments, performed in [5] showed that actually sin #13 7^ at 
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> 3a. The authors of [JJ find: 

sin 2 0i3 = 0.021 (0.025) ± 0.007 , 



(1) 



using the "old" ("new") fluxes of reactor v e in the analysis. Moreover, it was found in the 
same global analysis that cos 5 = — 1 (and sin #13 cos 5 = —0.14) is clearly favored by the 
data over cos 5 = +1 (and sin #13 cos 6 = +0.14 Q| ). 

The T2K and MINOS results will be tested in the upcoming reactor neutrino exper- 
iments Double Chooz [6], Daya Bay [7] and RENO [8j. If confirmed, they will have far 
reaching implications for the program of future research in neutrino physics. A relatively 
large value of #13 opens up the possibilities, in particular, i) for searching for CP violation 
effects in neutrino oscillations experiments with high intensity accelerator neutrino beams 
(like T2K, NOi^A, etc.); ii) for determining the sign of Amf^, and thus the type of neu- 
trino mass spectrum, which can be normally or invertedly ordered (see, e.g. [2]), in the long 
baseline neutrino oscillation experiments at accelerators (NOz^A, etc.), in the experiments 
studying the oscillations of atmospheric neutrinos (see, e.g. [9]), as well as in experiments 
with reactor antineutrinos [10]. A value of sin #13 ?Z 0.09 is a necessary condition for a 
successful "flavoured" leptogenesis when the CP violation required for the generation of 
the matter-antimatter asymmetry of the Universe is provided entirely by the Dirac CP 
violating phase in the neutrino mixing matrix |llj . 

In the present article we investigate the possibilities to obtain the relatively large value 
of #13, suggested by the T2K, MINOS and the current global neutrino oscillation data, in 
the case when the diagonalization of the neutrino mass matrix gives a negligible contri- 
bution to #13. Parametrically the latter means a contribution smaller than 0(A 2 ), where 
A = 0.23 is the Cabibbo angle expansion parameter. The diagonalization of the neutrino 
mass matrix leads to a negligible #13, for example, in models with L e — — L T flavour 
symmetry [12] and in most models accounting for bimaximal (BM) [13] or tri-bimaximal 
(TBM) |14] neutrino mixing. For our purposes, this is equivalent to setting the neutrino 
mass matrix contribution to #13 to zero, which we will do in our further analysis. The #13 
mixing angle can then be due to the contribution to the neutrino mixing coming from the 
diagonalization of the charged lepton mass matrix (see, e.g. [15] [16] [17] and the references 
quoted therein I). 

In particular, we will consider the possibility of a "large" #13 originating from the 
charged lepton mass matrix in the context of SU(5) supersymmetric unification. The 
motivation for such a choice are strong and well known. In our case, this choice will also 
allow to get a handle on the expected contribution to #13 and on the structure of the 
charged lepton Yukawa coupling matrix. This is because the SU(5) symmetry allows to 
relate the latter and the down quark one, on which further experimental information is 
available. 

This result can have important implications for the "flavoured" leptogenesis scenario of generation of 
the baryon asymmetry of the Universe [5]. 

2 For alternative possibilities leading to a relatively large #13, see [18] as well as, e.g. |19] . 
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In minimal SU(5) models, the down-quark and transposed charged lepton mass matrices 
are equal. On the other hand, depending on the (renormalizable or non-renormalizable) 
SU(5) operator(s) giving rise to a given Yukawa coupling matrix entry, the relation may 
involve Clebsch-Gordan (CG) coefficients. Such coefficients need to be invoked in order to 
fix the SU(5) predictions for the rrib/m T and rris/m^ mass ratios (unless specific values of 
tan/3, peculiar patterns of SUSY breaking parameters, or highly asymmetric mass textures 
are employed [20]). as shown by the authors of [21] 122] after having identified the possible 
GUT scale mass ratios in a bottom-up approach. All the possible CG coefficients in SU(5) 
and the Pati-Salam group entering such relations arising from dimension five and some 
dimension six operators are listed in |22| . Such coefficients not only can give rise to new 
mass ratios, but can also affect the prediction for #13 |22| . In this paper, we discuss possible 
combinations of CG coefficients, which can give a "large" $13 even if there is no contribution 
to #13 from the neutrino mass matrix at all. 

The paper is organised as follows. In Section 2 we describe the general set up, the 
determination of the standard neutrino parameters under the assumptions we make, the 
numerical procedure we used to determine the allowed CG coefficients and the predictions 
for #13. The results are illustrated in Section 3. In Section 4 we consider the special case 
in which the neutrino contribution to the PMNS matrix is of the BM or TBM form, in 
which interesting implications for the leptonic (Dirac) CP-violation phase can be drawn. 
In Section 5 we summarize. 



2 General Setup and Procedure 
2.1 The Setup 

We start with some definitions and a summary of the assumptions we make regarding the 
structure of the mixing matrices. Thereby we follow closely the discussion in [161 ESI ITT] . 
The PMNS lepton mixing matrix U is given by 

T T T diag T j 

U = U e U$, with U e , U u defined by % E y \ (2) 

m u = Ujm^Uv 

where toe and m v are the charged lepton and light neutrino mass matrices, respectively, 
and m E ias , m^ iag are diagonal and positive. The PMNS matrix can be parameterized in the 
standard way as U = -R23(023)-Ri3(0i3j £)-Ri2(0i2)Qi where Q is a diagonal phase matrix 
containing the two physical Majorana CP violation phases [231 124] . 

Barring correlations among the entries of tue, the hierarchy of the charged lepton 
masses translates into the possibility of diagonalizing uie perturbatively through subse- 
quent 2x2 unitary rotations in the following order: U e = U^U^U^, where Ufj is a unitary 
rotation in the ij block. The 13 rotation is negligible for our purposes and will be set to 
zero, thus Uf 3 = 1. 
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The assumption of zero neutrino mass matrix contribution to #13 can be rephrased as 
#13 = in the limit of U e = 1, or equivalently as Uj, = U\ije =1 = R2 3 (6 23 )Ri2(0'( 2 ) U P to 
phase matrices. Note that since U 23 does not contribute to #13, the assumption can be 
further rephrased as 



U = Ri 2 (Ot 2 )$R 23 (0 23 )R 12 (0 12 )Q 



C12C12 - s 12 c 23 sl 2 e^ 
-C12S12 - Sl2C23Cf 2 e i 

S12S23 



Sl2Cf 2 + Cl 2 C23S?2 ei 



h3S e 12 e^ 



Cl2C23q 2 e' v - Si 2 Sf 2 5236126" 
-C12S23 C 2 3 



(3) 



Here $ = Diag(l, e^, 1), where is a CP- violating phase [15] • The angles 9f 2 , 9±2, 623 
and the phase <j) in the expression ([3]) for U can be arranged to lie in the intervals [0, 7r/2] 
and [0, 27r], respectively. 

By commuting with the 623 one, the #i 2 -rotation has induced a non vanishing #13 given 

by 

sin 6f 2 tan #23 



sin 13 



U, 



1 + tan 2 6 



23 



sm 



2ae 
7 12 



smt/ 12 sm t/23. 



(4) 



Note that we have traded here 623 with #23 in the standard parameterization. Note also 
that the corrections to 623 and #12 induced by 6f 2 are given by 



tan #23 = cos 6*12 tan 9 23 

1 + cos #23tan6»f 2 e^ 



tan 6\ 2 = tan #12 



1 — tan #12 cos 623 tan 8f 2 e 



(5a) 
(5b) 



In the case of 3-v mixing under discussion, the magnitude of the CP violation effects 
in neutrino oscillations is determined [25] by the rephasing invariant associated with the 
Dirac CP violating phase 5: 



j C p = im{u: 1 u; 3 u e3 u„ 1 ] . 



(6) 



The rephasing invariant Jcp, as is well known, is a directly observable quantity. It is analo- 
gous to the rephasing invariant associated with the Dirac phase in the Cabibbo-Kobayashi- 
Maskawa quark mixing matrix, introduced in [26]. In the standard parametrisation of the 
PMNS matrix we find 



Jcp = - sin 26*23 sin 26>i2 cos #13 sin 26*13 sin 5 . 



(7) 



Using eqs. ([6]) and 



Jcp 



we get [T7] 
1 



sin 26^2 sin 2#i2 sin26>23 sin 623 sin< 



(8) 
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Operator Dimension 

4 1 
-3 

5 -1/2 

1 

±3/2 
-3 
9/2 
6 
9 

-18 



Table 1: Summary of possible SU(5) predictions for the coefficients ay, Numbers are 
taken from |22| . where also the corresponding operators are listed. 

Thus, to leading order in sin#f 2 we have sin 5 = — sin0. By comparing the real part of 
^ei ^3 m the two parameterisations we conclude that, at the same order, 

6 = -</>. (9) 
2.2 Relation between 0^ 2 and 9 13 in GUTs 

We would like now to study the possibility to generate a 6f 2 large enough to induce a #13 
in the range indicated by recent experiments in the context of an SU(5) Grand Unified 
Theory (GUT). The unification assumption is powerful because it allows to relate the 
charged lepton and down quark Yukawa matrices A^ and Aj> If all the Yukawa entries 
were generated by renormalizable operators and the MSSM Higgs fields were embedded in 
5 and 5 representations only, we would have = \®, leading to wrong predictions for 
the fermion mass ratios. In the general case one has instead A^ = cVjjA^. The coefficients 
a.ij depend on the operators from which the Yukawa entries arise. Such values can be 
constrained to belong to a finite set of rational numbers at the price of assuming that each 
Yukawa entry comes at least dominantly from a single renormalizable or non-renormalizable 
SU(5) operatoid. In this case, the possible values of the coefficients are listed in Tabled! 
see also [22]. 

The 9\ 2 angle is obtained from the diagonalization of the 12 block of the charged lepton 
Yukawa matrix after the 23 block has been diagonalized. Let us denote such 12 blocks in 
the charged lepton and down quark sectors (in the RL convention in which the Yukawa 

3 This could not be the case, for example, if SU(5) is embedded in SO(10) or a larger unified group. 
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interactions are written with the left-handed fields on the right) as 



A [12] 



a b' 
b c 



A 



[12] 



cea 



Bb 
7 C 



(10) 



In the following we will assume that the entries in eq. (|10p can be approximated with the 
corresponding entries of X E ' D , in which case the coefficients a, (3, f3', 7 are still bound to 
take one of the values in Table Q] (the rotation used to diagonalize the 23 sector can have a 
sizeable effect on the coefficient 7 and, if the charged lepton contribution to #23 from Ufa 
is sizeable, on the coefficient (3). 

We would like to determine the values of the coefficients a, /3, B' , 7 allowed by data, 
and in particular capable to account for the indication for a sizeable #13. Not all the values 
of the coefficients are allowed, in principle. The observables to be described are in fact 



'l3 j \Vus\i 



m. 



m,. 



m. 



(11) 



and for given a, (3, f3' , 7, the five experimental inputs above depend on the four real 



variables \b/c\, \b'/c\, \a/c\, lo, where the phase uj is defined by ac(bb')* 
explicit dependence is given by the following relations 



tan 6 



12 

V us \ 
m e 

m d 
m s 
m 



7c 

' 1 1 



1 



f3b 
7c 
b' 2 



a a 
7 c 
a bb' 



c 

(3(3' bb 1 



Bb* aa 
7c* 7c 
b 2\ b'* 



a 

+ — - 

c* c 



7 2 c 2 



1 - 
\b\ 2 



± A 

f3 2 \b\ 2 + f3' 2 \b'\ 2 



7 2 |c| 2 



h'\2 



^ = W 1 + 



2|c| 2 7 2 



U>\2 



e lu \acbb'\. The 

(12a) 
(12b) 
(12c) 
(12d) 
(12e) 



where A takes into account the possibility of a model-dependent contribution to | V us 

0.045. 



from 
The 



the up quark sector and is assumed to be in the range |A| < \J m u /m, 
experimental inputs used for the quantities on the LHS are listed in Table El The relations 
above are approximated and are accurate up to corrections of order A 4 , if \b/c\ < \t>' /c\ < A, 
\a/c\ < A 2 . 

Besides the general case in eq. (|10p . we will also consider the case in which a = and 
the symmetric case in which |A^| = IA21I and |Af|| — |A^|, as they arise in many models 
of fermion masses. Note that the symmetry condition implies b = ±6' and (3 = (3' . 
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Input Parameter 


Value 


Assumed error distribution 




(4.7362 - 4.7369) x 10" 3 [29] 


Uniform 




2.48 - 7.73 [22] 


Uniform 


m s /m d 


18.9 ±0.8 [3D] 


Gaussian 


\v us \ 


0.2252 ± 0.0009 [M] 


Used uniform in \V US \ ± A 


sin 6»i3 


0.089 - 0.190 (2a) [4] 


Approximated with a Gaussian 



Table 2: List of input parameters used in our analysis. 



2.3 Procedure 

Before we come to the results we briefly discuss the procedure we implemented. Since only 
the ratios \a/c\, \b/c\, and \b'/c\ enter when computing the experimental inputs, we have 
set in our forthcoming numerical analysis |c| = 1. We can also always perform a phase 
redefinition of the fields such that all the remaining coefficients are real and positive and 
the only physical phase is in a, so that a = exp(icj). 

For each possible combination of CG coefficients we diagonalized exactly both the mass 
matrices, using the expressions for the observables in eq. (Ilip in terms of a,b,b',u, of 
which the relations in eq. ()12[) are the expansion at NLO. Then we determined numerically 
a solution for these parameters such that all the experimental inputs are satisfied. We 
extracted the values for these inputs randomly following the distributions given in Table [2j 
We repeated this procedure until one solution is found. If, after a large number of attempts, 
no solution is found, we discard this combination of CG coefficients. For the viable CG 
coefficients we obtained by this procedure a distribution for #13, from which we computed 
the mean value and the standard deviation. To obtain sin #13 from sin 9f 2 we have assumed 
that #23 in the neutrino sector is maximal for simplicity. Given the uncertainties on the 
other input variables, this is a good approximation. 

Note that eq. (|12e|) fixes 7 to lie in the range of the observed m M /m s . Therefore we 
used this equation only to reduce the possible values of 7 to —3, 9/2 and 6, cf. Table [2j 
The GUT scale ratio m^/m s depends strongly on low energy SUSY threshold corrections 
[27] and in principle one can use them to push this ratio to more extreme values, but in 
simple SUSY breaking scenarios these are the only plausible values [22]. 

3 Results 

In the most general case it is easy to obtain values of the CG coefficients leading to a value 
of sin #13 compatible with the recent fits in [3J. In fact we find several hundred possible 
combinations. We therefore restrict ourselves in the following to some well motivated cases. 
We give a graphical summary of the results at the end of this section in Fig. [TJ 



7 



{a, 0,0', 7} 



{a, 6, 6', a;} 



sin 6»i3 



{-3,1,-3,-3} {0.0151,0.220,0.189,-2.81} 0.130 ±0.013 



Table 3: Possible CG coefficients with Yukawa couplings coming only from rcnormalizablc 
operators. We also show typical values for the entries of Ad, where c is normalised to one, 
and we give the prediction for sin #13 inlcuding its la standard deviation. 



{a, p,?, 7} 


{a, b, b',Lu} 


sin 6>i3 


{1,-3/2,-3/2,6} 

{1,-3/2,6,6} 

{-3/2,1,6,6} 

{6,1,6,6} 

{6,-3/2,6,6} 


{0.0899,0.246,0.679,0.145} 
{0.0286,0.212,0.153,-2.57} 
{0.0224,0.217,0.186,-2.34} 
{0.0155,0.281,0.259,0.278} 
{0.0134,0.247,0.184,-2.77} 


0.114 ±0.014 
0.103 ±0.008 
0.122 ±0.015 
0.175 ± 0.009 
0.137 ±0.014 



Table 4: Possible CG coefficients with Higgs fields in representations not larger than the 
adjoint. We also show typical values for the entries of Ad, where c is normalised to one, 
and we give the prediction for sin #13 inlcuding its la standard deviation. 

3.1 Results for Renormalizable Operators Only 

We start our discussion with the case in which the Yukawa couplings come only from 
renormalizable operators. This case is very restrictive as there are only two possible CG 
coefficients, which are = 1, if the Higgs sits in a 5 of SU(5), and otij = —3, if the Higgs 
sits in a 45 of SU(5) [28]. 

There is only one combination which is in agreement with the experimental data. It is 
shown in Table [3l where we give in addition typical values for the entries of Xd and the 
prediction for sin #13. 

3.2 Results without Representations larger than the Adjoint 

The next case we consider is the one in which the Yukawa couplings are generated by 
a dimension five operator, with all fields sitting in a representation not larger than the 
adjoint. This concerns also the messenger sector of a possible UV completion. Especially 
the Georgi-Jarlskog factor of —3 [28] is here not possible anymore. There are only three 
ciij left, which are 1, —3/2, and 6, giving five valid combinations as listed in TableHl where 
we give again typical values for the parameters and the predictions for sin #13, including 
its standard deviation. 

It is interesting to note that this possibility, as well as the last possibility, can be ruled 
out not only by a precise measurement of the leptonic mixing parameters, but also by a 
measurement of the SUSY spectrum. A CMSSM like spectrum with a positive fi parameter 
prefers a ratio m^/nig in the region of 4.5 — 6 |2 1 1. 122 j . ruling out the special case mentioned 
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{/M', 7} 


{b, b'} 


sin 6>i3 


{-1/2,-3/2,-3} 


{0.217,0.267} 


0.094 ± 0.003 


{-1/2,3/2,-3} 


{0.216,0.268} 


0.094 ± 0.003 


{-1/2,6,6} 


{0.251,0.240} 


0.164 ±0.013 


{1,-3,6} 


{0.212,0.273} 


0.094 ± 0.003 



Table 5: Possible Clebsch Gordan coefficients with a texture zero in the 11 element, a = 0. 
We also show typical values for the entries of A_d, where c is normalised to one, and we give 
the prediction for sin #13 inlcuding its la standard deviation. 

before. To get a small ratio m^jm s ~ 3, a spectrum more similar to an AMSB like scenario 
is preferred, in which the sign of the QCD part of the SUSY threshold corrections is flipped 
compared to the CMSSM with [i > 0. 

3.3 Results for a = 

The next scenario we discuss is a scenario, where we have a texture zero in the 11 element, 
a = 0. This can be motivated by having a flavon vacuum alignment, which has a zero 
in this position or having a Froggat-Nielsen mechanism at work, which puts there a zero 
or suppresses this element very strongly. For the CG coefficients we take all the possible 
values in Table [TJ In this case we end up with four possible combinations, which are listed 
in Table [5j Note that in this case there are no physical phases. 

3.4 Results for Symmetric Mass Matrices 

In the (anti-)symmetric case |Af^| = |A^| and \Xf 2 \ = (Afjj, which implies b = ±b' and 
P = (3', we find 26 possible combinations listed in Table [6l Such a mass matrix is generated, 
if the 12 and the 21 entries are coming from the same operator. Note that by choosing the 
unphysical phases appropriately we can always make b = b' . 

This case cannot be combined with any other case. If we restrict ourselves to certain 
operators or choose a = 0, no combination remains viable. 

4 Possible Implications for Dirac CP Violation in the Lepton 
Sector 

The results of the T2K and MINOS experiments [IJ [3] and of the global analysis of the 
neutrino oscillation data [4] have important implications for the Dirac CP violation in the 
lepton sector when the unitary mixing matrix Ul, originating from the diagonalisation of 
the neutrino mass matrix, is (up to CP violating diagonal phase matrices) of BM or TBM 
form, while that originating from the diagonalisation of the charged lepton mass matrix 
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{a,/3,7} 


{a,b,u} 


sin 13 


{-1/2,-3/2,-3} 

{-1/2,3/2,-3} 

{-3/2,-3,-3} 

{-3,-3/2,-3} 

{-3,3/2,-3} 

\0, —6, —6} 

{9,-3,-3} 

{-18,-3/2,-3} 

{-18,3/2,-3} 

{-18,-3,-3} 


{0.122,0.259,-0.183} 
{0.125,0.255,-0.0985} 
{0.115,0.233,-0.0736} 
{0.0205,0.284,0.243} 
{0.0143,0.268,0.201} 

•[U.UloD, U.ZUD, — o.Uoj 

{0.0142,0.212,-3.04} 
{0.0028,0.257,-0.294} 
{0.0033,0.255,-0.187} 
{0.0120,0.268,-0.0757} 


0.0903 ± 0.0008 
0.0903 ± 0.0008 
0.164 ±0.007 
0.098 ± 0.002 
0.098 ± 0.002 

u. lo» ± U.UUl 

0.144 ±0.003 
0.0901 ± 0.0008 
0.0900 ± 0.0007 

0.183 ±0.004 


{1,-3,9/2} 

{3/2,9/2,9/2} 

{6,-3,9/2} 

{9,-3,9/2} 

{-18,-3,9/2} 

{-18,9/2,9/2} 

{-18,6,9/2} 


{0.115,0.234,0.195} 
{0.093,0.186,0.107} 
{0.0254,0.289,0.108} 
{0.0155,0.275,0.137} 
{0.0057, 0.240, 2.97} 
{0.0117,0.209,-3.05} 
{0.0183,0.258,-3.10} 


0.105 ±0.009 
0.128 ±0.003 
0.1325 ±0.0009 
0.129 ±0.003 
0.107 ± 0.002 
0.149 ±0.003 
0.184 ±0.002 


{1,-3,6} 

{1,9/2,6} 

{3/2,9/2,6} 

{9/2,-3,6} 

{9/2,9,6} 

{6,-3,6} 

{9,-3,6} 

{-18,9/2,6} 

{-18,6,6} 


{0.127,0.258,-0.114} 
{0.0888,0.181,0.200} 
{0.111,0.225,0.177} 
{0.0200,0.280,0.105} 
{0.0918,0.179,0.060} 
{0.0182,0.280,0.249} 
{0.0108,0.263,0.286} 
{0.0094,0.221,-2.96} 
{0.0133,0.211,-3.08} 


0.0903 ± 0.0008 
0.097 ± 0.002 

0.11 ±0.01 
0.100 ±0.001 
0.183 ±0.002 
0.098 ± 0.002 
0.094 ± 0.003 
0.116 ±0.002 
0.143 ±0.003 



Table 6: Possible Clebsch-Gordan coefficients with a symmetric mass matrix and the 
resulting prediction for sin #13 . 
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rm n 



0.12 0.14 0.16 



0.18 



sin 8\ 



Figure 1: Graphical presentation of the results obtained. The vertical straight (dashed) 
lines denote the la (2a) allowed ranges of sin #13 taken from [3] (see also Table [7]). The 
yellow, green, blue, red bins correspond to the results from Table [3J HI [SJ [51 



can be approximated as U e = ^12(^12) ( see eq. ([3])). We recall that if U$ coincides with 
the BM mixing matrix one has: 



sm 



'12 



1 

2 ' 



sm 



'23 



1 

2 ' 



sm 



'13 



. 



(13) 



rt 



sm 



7 13 



. 



(14) 



In the case of tri-bimaximal mixing form of Ul, 

sin 2 0fc = I, sin 2 0^ = 1 , 

While in both cases 2 3 coincides, or is close to, the experimentally determined best fit 
value of the atmospheric neutrino mixing angle 623, in the case of bimaximal mixing form 
of ul a relatively large correction to d\ 2 is needed to get a value of the solar neutrino 
mixing angle #12 compatible with that determined from the data. And in both cases a 
non-zero #13, having a value in the range sin 2 #13 = 0.021 (0.025) ± 0.007, see also eq. ([II), 
suggested by the current data, has to be generated. 

If Uu has bimaximal form, the angle Q\ 2 = ?r/4 is corrected by the charged lepton 
mixing as follows: 



sm f 12 — — 



1 1 + cos (j) sin 20f 2 ~ \ sm2 #12 



1 - \ sin 2 0f 2 



(BM) 



(15) 



where we have neglected the possible contributions from the charged lepton mixing angles 
<9f 3 and 0f 3 0. To leading order in sin(9f 2 we get [El Q31 El [32] : 



1 1 



1 



sin 0i2 — - H — p cos (f) sin 0f 2 — 77 + cos ^ s i n #13 > (BM) , 
2 y 2 2 



(16) 



'These contributions are given in [151 117| . 
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sin 13 



sin 2 i2 sin 2 23 



Best fit 0.145 0.306 0.42 

la 0.114 - 0.167 0.291 - 0.324 0.39 - 0.50 

2a 0.089 - 0.190 0.275 - 0.342 0.36 - 0.60 

3a 0.032 - 0.210 0.259 - 0.359 0.34 - 0.64 



Table 7: Results of the global fit of the PMNS mixing angles taken from [4] and used in 
our analysis. The results quoted were obtained using the "old" reactor v e fluxes (see [4] for 
details). 



where we have used the relation cos (f> = cos 5 and the fact that in the approximation 
employed we have: 

sin 0i3 ~ -=sm9f 2 > 0. (17) 
V2 

The sign of the second term in the r.h.s. of the equations in (|16|) is important H in view of 
the fact that, according to [I], the global neutrino oscillation data favors a negative value 
of cos 5 sin 0i3 over the positive value. We note that such a sign is in agreement with the 
one needed to reduce the bimaximal prediction for 0i2 down to the experimentally allowed 
range. 

In the case of the TBM form of ut we find under the same assumptions: 

11 + cos d> sin 20? 



sin 2 0i2^-— 2f)e 12 , (TBM). (18) 

o 1 — ± sm z 0f 2 

To leading order in sin0f 2 we get: 



,2„ 1 , ^ 



sin 2 0i 2 ~ - + -1- cos S sin 13 , (TBM) , (19) 

where we have used again cos(/> = cos 5 and eq. (|17p . which is valid also in this case. 
The angle 23 = vr/4 gets the same correction in both cases of BM and TBM lit: 



1 cos 2 0? 
2~T 



sin 2 23 ^ —I 2 . (20) 



7 12 

As it follows from the above expression, the leading correction is of order sin 2 0^ 2 [IS 



We note that in [151 1171 [32] a somewhat different parametrisation of the PMNS matrix was used 
(namely, U = U\U V and U = R12 (^12)^^23 (^23 )Ri2 (Q12 )Q) and instead of the first relation in (fTo) . the 
relation sin 2 612 — 1/2 — (eos</> sin&i 2 )/\/2 was obtained. It is not difficult to show that in this case we 
have cos S = — cos(f> and thus one arrives at the same result for the relation between sin 2 812, 613 and S. 
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0.00 0.05 0.10 0.15 0.20 0.00 0.05 0.10 0.15 0.20 

sin #i3 sin #13 



Figure 2: The cosine of the Dirac CP phase S as a function of sin $13 in the cases of 
tri-bimaximal (left panel) and bimaximal mixing (right panel) arising from the diagonali- 
sation of the neutrino mass matrix, eqs. (TT5)) and (TT5)) . The green, yellow, orange regions 
correspond to the 1, 2, 3er allowed ranges of sin 2 #12. In the bimaximal mixing case only 
the negative values of cos S are compatible with the data on sin 2 # 12 . The vertical straight, 
dashed, dotted lines denote the 1, 2, 3a allowed ranges of sin# 13 . The values of sin 2 # 12 and 
sm6> 13 are taken from [3] (see text for details). 



The values of sin 2 6*12, sin 2 #23 and sin 2 #13, obtained in the global data analysis in 0] 
are shown for convenience in Table In Fig. [2] we present graphically the constraints on 
cos 5 implied by the data on sin 2 812 and sin #13 in both cases of BM and TBM ut, using 
the relations (I15( [T8j) , which are exact in 8\ 2 • 

A few comments are in order. It follows from our analysis that in the case of bimaximal 
ut, the current 3a (2a) experimentally allowed range for sin 2 8u requires that sin #13 <; 
0.14 (0.16), the minimal value being very close to the best fit value found in [3]. If future 
data on #13 will show (taking into account all relevant uncertainties) that, e.g. sin #13 ^ 0.10, 
the simple case under discussion of ut having a BM mixing form and U e = ^12(^12) wn l 
be ruled out. Further, using the 3a (2a) allowed ranges of both sin 2 812 and sin #13 we 
find that cos 6 is constrained to lie in the interval: —1 < cos 5 ^ —0.60 (—0.79). Thus, 
cos 5 = 0, or 5 = n/2, and therefore maximal CP violation in neutrino oscillations, is ruled 
out in the scheme we are considering. 

We get very different constraints on cos 5 in the case of the TBM form of ut- If we 
use the 3a allowed ranges of sin 2 812 and sin #13 in the analysis, all possible values of 
cos 5 are allowed. The 2a intervals of allowed values of sin 2 812 and sin #13 require that 
—0.66 ^ cos 5 ^ 0.09. Thus, maximal CP violation, 5 = it/2, is allowed. 

Using the constraints on 5 derived above, we have obtained also predictions for the 
magnitude of the rephasing invariant Jcp- The latter has the following form to leading 



13 



0.02 



y CP o.oo 



-0.04 




0.04- 



0.02 - 



-0.04 - 



0.10 
sin #13 



Figure 3: The rephasing invariant Jcp as a function of sin #13 in the cases of tri-bimaximal 
(left) and bi-maximal (right) mixing from the neutrino sector. The green, yellow, orange 
regions correspond to the 1, 2, 3tr allowed ranges of sin 2 6*12. The vertical straight, dashed, 
dotted lines show the 1, 2, 3cr allowed ranges of sin #13. The 1, 2, 3a allowed ranges of 
sin 2 9 V2 and sin # 13 are the same as in Fig. [2] 



order in sin in the cases of BM and TBM Ut: 



BM 
CP 



— sin 5 sin #13 , 



tTBM 

j c~p 



3V2 



sin 5 sin 9 
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(21) 
(22) 



The results of this analysis are presented in Fig. [3l As could be expected on the basis 
of the results obtained for cos S, the existing data on sin 2 #12 and sin #13 imply that in the 
case of BM form of lit, Jqp should lie in the interval | Jcp| ^ 0.037(0.031), using the 3a 
(2a) ranges, the value Jcp = being allowed. We get similar results for the maximal 
possible value of |Jcp| if Ut exhibits TBM mixing form. However, if we use the 2a allowed 
ranges of sin 2 9\2 and sin #13, we find that |Jcp| has to be non-zero and not smaller than 
0.014: 0.014 < I Jcp I < 0.037. Values of |J C p| k, 0.01 are potentially accessible by the 
accelerator experiments T2K, NOi^A, etc. planning to search for CP violation effects in 
neutrino oscillations (see, e.g. |33j). 



5 Summary and Conclusions 

In the present article we have considered the possibility that a sizeable value of the neutrino 
mixing angle #13, compatible with the recent indications from the T2K, MINOS and the 
global neutrino oscillation data, arises primarily from the contribution of the charged lepton 
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sector U e to the lepton mixing, U = U e Ul, where U is the PMNS matrix and Ul arises 
from the diagonalisation of the neutrino mass matrix. This scenario should necessarily 
be taking place in the context of a wide variety of models predicting a relatively small 
contribution to #13 from the neutrino sector, say #13 < 0(A 2 ), where A is the Cabibbo 
expansion parameter A ~ 0.23. The analysis has been performed in the context of SU(5) 
supersymmetric models, which are independently well motivated and allow the charged 
lepton contribution to #13 to be related to the mass and mixing observables in the quark 
sector. 

The naive relation mg = between the charged lepton and down quark mass matrices 
is not compatible with the indication of a relatively large #13. This is not surprising, as 
the quoted relation is not compatible with the well known values of the charged fermion 
mass ratios either and it has to be corrected, possibly by Clebsch-Gordan (CG) factors. 
This is particularly true in the case (relevant for us) of the charged fermions belonging 
to the first two families, whose small masses can arise from non-renormalizable operators 
that can well involve SU(5) breaking sources. Under motivated hypotheses for the values 
of such CG coefficients, we found all the combinations that are experimentally viable. The 
main results are summarized in Tables El [H (U and [HJ In order to describe the present 
data, one needs either a non- vanishing 11 entry in the down quark and charged lepton mass 
matrices, or CG coefficients from non-renormalizable or from SU(5) breaking operators in 
representations larger than the adjoint. 

We also studied the specific case in which the contribution from the neutrino sector to 
the lepton mixing, Ul in eq. ([2]), is, up to the phase matrices $ and Q, see, eq. ([3]), of the 
bimaximal (BM) or tri-bimaximal (TBM) form, while that from the charged lepton sector 
is assumed to be of the form U e = R^i^^)- Under these assumptions, $ = diag(l, e 1 ^, 1), <p 
being a CP violating phase, and Q contains the two Majorana CP violating phases. In this 
case, it is possible to draw important conclusions on the leptonic (Dirac) CP-violating phase 
5 ?» —(f) and on the magnitude of the CP violation effects in neutrino oscillations. The latter, 
as is well known, is determined by the rephasing invariant Jcp- Our results are summarized 
in Figs. [21 [3j We find that in the case of bimaximal Uv, the current 3a (2a) experimentally 
allowed range for sin 2 9\2 requires that sin #13 ^ 0.14 (0.16). If future data on #13 will 
show (taking into account all relevant uncertainties) that sin #13 has a smaller value, e.g. 
sin #13 ^ 0.10, the simple case of U v having a BM mixing form and U e = Ri2(0f 2 ) will be 
ruled out. Further, using the 3a (2a) allowed ranges of both sin 2 #12 and sin #13 we find 
that cos S is constrained to lie in the interval: —1 < cos 5 ^ —0.60 (—0.79). Thus, cos 5 = 0, 
or 5 = it/2, 37r/2, and therefore maximal CP violation in neutrino oscillations, is essentially 
ruled out in the scheme we have considered. If it will be confirmed experimentally that 
cos 5 sin #13 < 0, that would imply (in the standard parameterization of the PMNS matrix, 
in which sin6»i 3 > 0) that 127° (142°) < 5 < 233° (218°). Correspondingly, J CP should lie 
in the 3a (2a) interval |Jcp| < 0.037 (0.031), the value Jcp = being allowed. However, 
if it will be experimentally established that cos 5 sin #13 > 0, that would rule out the case 
of U e = #12(012) and Ul having a BM mixing form. 
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In the case of the TBM form of Ul and U e = i?i2(# 12 ), all possible values of cos (5 are 
allowed if one uses the 3a allowed ranges of sin 2 9±2 and sin #13 in the analysis. The 2a 
intervals of allowed values of sin 2 9\2 and sin #13 require that —0.66 cos<5 0.09, or 
85° ^5 5 Si 131°. Thus, maximal CP violation, 5 = vr/2, in both cases is allowed. We 
find using 2a allowed ranges of sin 2 612 and sin #13, that |Jcp| has to be non-zero and not 
smaller than 0.014: 0.014 < | Jcp| ^$ 0.037. This interval falls in the range of the potential 
sensitivity of the T2K and the future neutrino neutrino oscillation experiments (NOi^A, 
etc.), designed to search for CP violation effects in neutrino oscillations. 
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Note Added. During the finalising stage of the work on the present article, Ref. [33] 
appeared, where the authors follow a similar approach in generating a relatively large #13. 
Our work goes beyond the analysis performed in |34| in a number of aspects. In [34] only 
the case a = is discussed, while we consider the more general case of a ^ 0. The authors 
of [M] use only leading order terms (without error estimates) in the expressions for the 
observables in terms of the parameters of the down quark and charged lepton mass matrices 
(i.e., in eqs. (12a) - (12e)) and they do not utilise the rather precisely known value of the 
ratio m s /md (see Table [2]). This rules out some of the possibilities considered in [34] to 
be viable. Finally, the authors of |34| draw only qualitatively conclusions about the value 
Dirac CP violating phase 5 in the case of (tri-)bimaximal mixing arising from the neutrino 
sector. We perform a quantitative analysis and determine the allowed ranges (at 1,2 and 
3 a) of values of cos 5 and of the rephasing invariant Jcp in the two cases (Figs. [2]and[3j). 
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